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ABSTRACT 

Context. Acoustic waves are one of the primary suspects besides magnetic fields for the chromospheric heating process to tempera- 
tures above radiative equilibrium (RE). 

Aims. We derived the mechanical wave energy as seen in line-core velocities to obtain a measure of mechanical energy flux with 
height for a comparison with the energy requirements in a semi-empirical atmosphere model, the Harvard-Smithsonian reference 
atmosphere (HSRA). 

Methods. We analyzed a 1-hour time series and a large-area map of Call H spectra on the traces of propagating waves. We analyzed 
the velocity statistics of several spectral lines in the wing of Ca II H, and the line-core velocity of Ca II H. We converted the velocity 
amplitudes into volume (oc pv 2 ) and mass energy densities (oc v 2 ). For comparison, we used the increase of internal energy (oc RpAT) 
necessary to lift a RE atmosphere to the HSRA temperature stratification. 

Results. We find that the velocity amplitude grows in agreement with linear wave theory and thus slower with height than predicted 
from energy conservation. The mechanical energy of the waves above around z ~ 500 km is insufficient to maintain on a long-term 
average the chromospheric temperature rise in the semi-empirical HSRA model. The intensity variations of the Ca line core (z ~ 1000 
km) can, however, be traced back to the velocity variations of the lowermost forming spectral line considered (z ~ 250 km). 
Conclusions. The chromospheric intensity, and hence, (radiation) temperature variations are seen to be induced by passing waves 
originating in the photosphere. The wave energy is found to be insufficient to maintain the temperature stratification of the semi- 
empirical HSRA model above 500 km. We will in a following paper of this series investigate the energy contained in the intensity 
variations to see if the semi-empirical model is appropriate for the spectra. 
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1. Introduction 

The chromospheric and coronal temperature rise above the pho- 
tospheric temperature requires a heating mechanism in addition 
to radiative energy transport in the outer atmospheres of stars. 
The existence of chromospheres and coronae on the Sun and 
other stars is a well established fact, but still no generally ac- 
cepted picture of the energy transport process responsible for 
thei r formation could be derived, especially for the corona (see 
e.g. iNarain & Ulmschneiderll 1996b . 

The scenarios suggested for the chromospheric heating pro- 
cess can be divided in three main categories: 1. field-free heat- 
ing by the propagation of (acoustic) waves in a gas plasma, 2. 
indirect magnetic heating in which magnetic field lines act as 
"catalyst" without being destroyed in the process, and finally, 3. 
direct magnetic heating in which the reconnection of magnetic 
field lines leads to a conversion of magnetic to thermal energy, 
thereby removing the field lines. 

In the first case, acoustic waves are generated by mate- 
rial (connective) mot ions as purely mechanical phenomenon 
(Musi elak et al.|[l994l) . The change of the atmospheric proper- 
ties in the vertical direction in a gravitationally stratified medium 
leads to a "steepening" of vertically propagating waves produc- 
ing a shock front. The energy of the propagating wave is dis- 
sipated in the ambient medium in and near the shock front, 



and hence, mechanical energy can be transported over some 
distance before conversion to thermal energy. Because of the 
stochastic nature of the excited acoustic wa ves, the energy trans- 
port process would not be continuous (e.g. Ulmschneider 1 97 ll: 



luimschneide r et al.ll 19781: ICarlsson & Ste in 1997). 
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For the second case of indirect magnetic heating requiring 
the presence of magnetic fields as a catalyst, two different sce- 
narios are possible in quiet Sun regions due to the specific topol- 
ogy o f the magnetic fields: the (more) vertical strong network 
fields (Solankif ri993l) and the (more) inclined weaker internet- 
work (IN) fields dLites et aDl 19961: iKhomenko etaDl2003l 120051: 
lOrozco Suarez et aD l2007|) . The network contains small-scale 
flux concentrations (some hundred kms) in the form of isolated 
flux tubes or flux sheets. These concentrated magnetic fields re- 
side in the intergranular lanes, because the interior of granules is 
swept free of magnetic fields by the directed convective flows as 
long as the ratio of mechanical to magnetic energy, /3 ~ p gas /B 2 , 
is above 1 . The spatial extension of the flux concentrations thus 
is strongly restricted in the photosphere by the converging gran- 
ular flows, but in the chromosphere the fields can expand be- 
cause of the reduced outside g as pressure a nd the lack of directed 
mass flows (iHammerl fl987l: iKeller et all 11990b ISolankil 119931: 
ISteiner et al.ll 19981) . This leads to a specific topology with small- 
scale, mostly vertical fields in the photosphere that get more in- 
clined in the chromosphere and fill the whole atmospheric vol- 
ume at some height. In t he layer of field expansion, a nearly hor- 
izontal canopy exists dJonesI 1 19851: ISchaffenbe rger et all 120061: 




Fig. 1. Normalized intensity contribution function of Call H in a 
LTE calculation. Red: center of contribution f unction. The green 
line was derived from phase differences dBeck et al.l i2008). 



iLites et all [2008). In this topology, two possibilities exist to 
transfer energy from the photosphere upwards in or near the net- 
work. The first option is that acoustic waves are generated out- 
side of the magnetic fields where some of them will propagate 
upwards as before. The deposit of their mechanical energy is 
then not triggered by the steepening in the gravitationally strat- 
ified medium, but the contact of the propagating wave with (or 
transition into) the inclined m agnetic field lines (IChitre & Davilal 
1 1 99 ll; iDavila & Ch itre 1991). The second option is the propaga- 
tion of waves inside the magnetic fields themselves. These waves 
can come in different modes (fast and slow magneto-acoustic 
waves, Alfven waves; see, e.g., Narain & Ulmschneider 1996). 
They propagate along the field lines which with the topology 
of solar quiet Sun fields would yield automatically a predomi- 
nantly vertically propagation. The deposit of energy then would 
be governed by the vertical variation of atmospheric properties 
inside the field concentrations or mode conversion processes. 
The source of the magnetic-field-related waves would again be 
the convective motions of the surrounding granulation or p-mode 
oscillations that constantl y buffet the magnetic field concentra- 
tions in the photosphere (Hind man & Jainll2008l) . The indirect 
magnetic heating by waves inside the fields would tend to affect 
a larger volume in the chromosphere than in the photosphere, as 
the magnetic field lines directly connect a large chromospheric 
volume with its small photospheric counterpart. 

For the third case of direct magnetic heating by re- 
conne ction, again two options are possible. Recent observa- 
tions dOrozcoSuarez et al.ll2Q07b ILites et al.|[2Q08b ISteiner et ail 
2008]) have shown that a large amount of weak (B < 
0.5 kG) horizontal magnetic flux is present inside of gran- 
ules that seems to be swept upwards with the granular up- 
flows. The existence of this magnetic flux was already sug- 
gested to explain the large differences of average magnetic 
flux between measurements using the Zeeman effect (sensi- 
tive to cancella tion of sub-resol ution opposite-polarity fields, 
around 20 G; Khomenkoet ID [2005) and the Hanle effect 
(60-100 G, insensitive t o cancellation; iBianda et al.l Q 999; 
iTruiillo Bueno et al.|[20Q4b . iMartfnez Gonzalez et al.1 (120081) in- 
voked turbulent stochastically oriented fields to explain the small 
amount of center-to-limb variation of li near and c i rcular po- 
larization signal amplitudes in IN areas. ILites et all (120081) re- 
cently increased the Zeeman-based estimate of the flux amount 
to around 60 G, using seeing-free data from the Hinode satel- 
lite. The interaction of the weak horizontal flux with the mag- 
netic canopy of the network or an increasing amount of recon- 
nection between the horizontal fields with time could lead to 
an energy deposit in chromospheric layers, where the magnetic 
field energy would have been produced previously in the pho- 
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Fig. 2. Average Call H spectrum with spectral lines 1 to 9 
marked by vertical lines and intensity ranges OW, MW, H-index, 
and IW marked by pluses (+). 

Table 1. Line parameters. 



Nr. 


Ion 


A [ 


Exct. 


log(gf) 


Transition 






[nm] 


pot. [eV] 








1 


CrI 


396.36900 


2.544 


0.66 


5G 6.0- 


5H 7.0 


2 


Til 


396.42700 


0.021 


-1.18 


3F 3.0- 


3F 4.0 


3 


Fel 


396.45530 


2.839 


-1.55 


3P 1.0- 


5P 2.0 


4 


Fel 


396.55088 


3.241 


-1.778 


5D 3.0- 


5D 4.0 




Fel 


396.59152 


2.425 


-3.550 


3P 1.0- 


5D 1.0 


5 


Fel 


396.60617 


1.605 


-1.64 


3F 2.0- 


3D 3.0 




Fell 


396.64260 


1.724 


-5.525 


4P 0.5- 


6D 0.5 


6 


Fel 


396.66300 


3.212 


-0.360 


5D 4.0- 


5F5.0 




Fel 


396.68100 


3.301 


-1.949 


3D 2.0- 


3G3.0 


7 


Fel 


396.74206 


3.295 


-0.41 


3H 4.0- 


1G3.0 




Fel 


396.79614 


3.235 


-0.87 


5D 3.0- 


3G 4.0 


8 


Call 


396.84900 


0.000 


-0.166 


2S 0.5- 


2P 0.5 


9 


Fel 


396.93000 


1.485 


-0.42 


3F 4.0- 


3F3.0 




Fel 


630.15012 


3.654 


-0.75 


5P 2.0- 


5D 2.0 




Fel 


630.24936 


3.686 


-1.236 


5P 1.0- 


5D 0.0 



1: wavelengths partly courtesy of J.Bruls, KIS (per sonal note), 
partly collected from the NIST Atomic S pectr a Database dMartin et al.l 
1998; Reader et al. 2002) and Nave^GD G223) 



tosphere or below (cf. Ilshikawa & T suneta 2009). The second 
option of magnetic heating would be the reconnection of net- 
work fields triggered by the constant buffeting and braiding of 
field li nes due to the stochastic motion of the photosph eric foot 
points dNarain & Ulmschneidedll996l:lPeter et al.ll2004l) . 

Recently, iFontenla et al.l (|2008) have suggested the Farley - 
Buneman instability as another possible source of chromo- 
spheric heating where cross-field motions of ionized material 
lead to an energy deposit in a small height layer. The horizon- 
tal magnetic fields again only act as a catalyst in the process, but 
the motion of the ionized plasma is thought to come mainly from 
convective motion rather than waves. 

From the observational point of view, any heating process 
related to traveling waves is better accessible for detection than 
for example the direct magnetic heating processes (granted that 
they exist). The best studied chromospheric spectral lines, for 
which the importance of heating by waves, or at least, the influ- 
ence of waves has been definitely established, are Call H and 
K at a round 397 nm and 393 nm, respectively dLinsky & Avrettl 
Il970h . The Call H line allows for example to trace the signa- 
ture of waves at different atmospheric heights due to its broad 
line wings where only a region of around 100-200 km vertical 
extent contributes to the intensity at each wavelength (cf. Fig.Q] 



or Fig. 5 of iRezaei et al.ll2008l) . ICarlsson & Steinl (1 1997b have 
demonstrated that the (transient) bright grains seen in the core 
of H and K are due to acoustic waves forming shocks. It is 
at present strongly debated if (high-frequency) acoustic waves 
can supply enough ener gy to maintain a chromosphere (cf. 
iFossum & Carlssonll2005L and the vivid discussion spawned by 
it). The energy of the acoustic waves is, however, usually derived 
indirectly from measurements of intensity variations observed in 
(br oad-band) i maging . 

iBeck et all ([2008) showed that the intensity of the emission 
peaks of Ca II H has a clear phase relation to the intensity in the 
line wing for a large range of oscillation frequencies, i.e. bright- 
enings i n the core are preceded by the same in the line wing 
(cf. also lLiu|[l974l) . They also showed that the Ca spectra of net- 
work locations show indications of transient events like the field- 
fre(fl IN regions, in addition to another contribution of steady 
emission. Thus, both magnetic network and field-free IN are sig- 
nificantly influenced by the presence of propagating waves. 

In this contribution, we want to investigate the energy of 
these propagating waves by a detailed analysis of spectra of Ca II 
H. The effect of the weak IN magnetic fiel ds on the emi s sion o f 
Call H was estimated to be negligible by IRezaei et al.1 (2007). 
As the weak (or no) fields cover the largest part of the solar sur- 
face in the quiet Sun, a limitation to wave heating with neglect 
of direct magnetic heating processes seems reasonable. Note that 
the indirect magnetic heating processes are still included if they 
rely on propagating waves that have a LOS component at disc 
centre. 

After the description of the data sets used (Sect. [2]), we de- 
termine the area fraction with signatures of magnetic fields in the 
chromosphere in Sect. to check whether we can ignore field- 
related effects on the velocity statistics. We then investigate the 
relation between the line-core intensity and velocity of a set of 
spectral lines in Sect. |U to exclude a convective origin for the 
patterns. We study the statistical properties of the LOS velocities 
and the relations between the spectral lines of different forma- 
tion height and compare the mechanical energy contained in the 
waves with requirements of a semi-empirical atmosphere model 
in Sect. The findings are summarized and discussed in Sect. [6] 
In a following paper of this series, we will investigate the energy 
in the intensity variations of the Ca II H line for comparison with 
the (mechanical) energy contained in the LOS velocities. 

2. Observations, line selection and overview maps 

We used two observations taken in the morning of 24/07/2006. 
The first is a time series of around 1 hour duration (see Beck 
et al. 2008 for details, around 30000 spectra); the second is a 
large-area map of 75" x ~ 70" extent (also around 30000 spec- 
tra). The data sets were obtained with the POlarimetric LIttrow 
Spectrograph (POLIS, IBeck et al.1 l2005bl) . In both cases, the 
Stokes vector around 630 nm was obtained together with inten- 
sity profiles of Call H. The slit width was 075, as was the step 
width of scanning. Spatial sampling along the slit was 0729. The 
integration time was 3.3 sec per scan step for the time-series, giv- 
ing a cadence of a r ound 2 1 sec (repeated small maps of 4 steps). 
As in IBeck et aD (120081) . we only used the last step with co- 
spatial spectra in 630 nm and Ca II H. For the large-area map, the 
integration time per step was twice as large (6.6 sec). The data 
were reduced by the res pective set of routi nes for flat field and 
polarimetric correction (iBeck et al.ll2005bllab . The noise level in 




Fig. 3. Overview over the 1-hour time series. Left column, top 
to bottom: continuum intensity at 630 nm (I c ), intensity of outer 
line wing of Ca II H (OW), middle wing (MW), inner wing (IW), 
H-index. Right column, top to bottom: integrated unsigned cir- 
cular polarization, LOS velocity of line 1, 5, 9 (blends in line 
wing), 8 (Call H line core). Black/white correspond to Doppler 
shifts of ±1 kms -1 for line 1, to +1.5 kms -1 for lines 5 and 9, and 
to =f5 kms -1 for line 8. For the locations marked red no line-core 
velocity of Call H could be defined; the red contour lines mark 
these positions in the Stokes V map. The white contour lines in 
v8 mark locations with strong V signal. The grid lines have a 
spacing of 10" (y-axis) and 7 min (x-axis), respectively. 



the polarization signal in continuum windows was around 1 -10 
of the continuum intensity for the time-series and 8-10" 4 for the 
large-area map. Using (pseudc0-)continuum intensity maps, the 
two wavelength ranges were aligned to be co- spatial. 

In the Ca II H data, we identified several of the photospheric 
spectral line blends in the line wing. We used a subset of 9 spec- 
tral lines with a sufficient line depth in our study (cf. Fig. [2]), 
together with the two Fel lines at 630.15 nm and 630.25 nm in 
the second channel. Table [T] lists all lines with line transition and 



1 "Field-free" is used here in the sense of "no detectable circular po- 
larization signal above noise level". 



2 For Ca II H a spectral region in the line wing without a photospheric 
blend was used. 





Fig. 4. Same as Fig.[3]for the large-area scan. The grid lines have 
a spacing of 10" in both axes. The blue rectangle marks an area 
of identical structure in I c and vl. 

adopted solar rest wavelengths. After deriving the velocity statis- 
tics of all of them, it became clear, however, that the lines 1 to 4 
(and 630.15 nm and 630.25 nm) had nearly identical properties, 
as had the lines 5 to 7. This can be related directly to their line 
depth, and hence, their similar formation height (cf. Figs. [T] and 
[2 for lines 1 to 4 and 5 to 7). We then decided to used only lines 
1, 5, 8, and 9 in the later evaluation. We attributed a "line-core" 
velocity to the displacement of the absorption core of Ca II H, 
but caution that its interpretation as a velocity is less secure than 
for the other spectral lines. It also turned out that the statistics of 
both time-series and large-area map were as good as identical, 
so we usually used quantities averaged over both data sets if not 
noted otherwise. 

Figures [3] and HI show overview maps of the observations. 
The right column shows the integrated unsigned Stokes V signal 
(at top) to mark locations with photospheric magnetic fields, and 
LOS velocities for the four chosen spectral lines. No reasonable 
velocities (v » 10 kms -1 ) could be defined for the Call H line 
core (line 8) on some locations because of the absence of a well- 
defined absorption core. These locations tend to be co-spatial 
to or close-by to strong photospheric magnetic fields and often 
show a very strong H2 V emission peak with a plateau of constant 



intensity to the red of it, making a definition of velocity arbitrary. 
The maximum polarization degree 

/? = maxVG 2 + t/ 2 + V 2 , (1) 

can be used to estimate the presence and amount of photospheric 
magnetic flux. The fraction of pixels with a polarization degree 
p above 1 % of I c that indicate concentrated magnetic fields, and 
which correspond also to the bright patches in, e.g. , the to p right 
panel of Fig.H is around 10 % (see lBeck & Rezaeil (I2009L A&A, 
in press) for a discussion of the dependence of p on the spatial 
resolution). 

3. Influence of chromospheric magnetic field 
topology on Ca II H spectra 

3.1. Indirect magnetic field signature in intensity 

The direct measurement of the chromospheric magnetic fields in 
the lower chromosphere has not been possible up to now for a 
number of reasons. The spectral lines like Call H that form in 
a suited height range of the solar atmosphere are very deep and 
broad, leading to a low light level in observations of high spa- 
tial resolution. The magnetic fields are much weaker than in the 
photosphere; they induce only a weak polarization signal via the 
Zeeman effect. A detection of these weak polarization signals 
requires a high signal-to-noise ratio that is only possible with 
long integration times that, however, strongly reduce the spatial 
resolution again. Even if POLIS was designed for polarimetry 
in 630 nm and Call H, the latter could only be done for so- 
lar features with str ong magnetic fields like pores or sunspots 
(see lBeck et al.ll2QQ5bb . Lacking the direct detection of chromo- 
spheric magnetic fields, often the bright structures seen in imag- 
ing observations of chromospheric lines are taken as tracers of 
the magnetic field lines. 

The map of the wavelength-integrated H-index in Fig. @l 
shows the strongest brightenings at locations with photospheric 
polarization signal. If one, however, reduces the averaged wave- 
length range significantly down to a few wavelength points, the 
intensity patterns change drastically. Figure \5\ shows the FOV 
of the large-area map as seen in the very line core of Call H 
(middle left). We averaged the spectra over a fixed wavelength 
range of ±3 pm around the line core at 396.849 nm. This in- 
tensity map shows several thin and elongated bright fibrils that 
originate from the patches with strong photospheric polarization 
signal. An identification of such fibrils with magnetic field lines 
is, however, not straightforward, if one compares with the line- 
core velocity map (top left). The intensity fibrils usually have a 
counterpart in the LOS velocity map, implying that the apparent 
increase of intensity in the map taken at a fixed wavelength is 
only due to a Doppler shift of the absorption core. We marked 
two prominent example of intensity fibrils with green parallel- 
ograms. The intensity fibrils have a counterpart in the velocity 
map with exactly the same morphology. This relation is not valid 
only for a few selected case, but throughout the whole FOV, as 
visualized by the yellow contours that mark strong blueshifts of 
the Ca line core. The shape of these velocity contours inside the 
green rectangle is complementary to the bright fibrils in the line- 
core intensity. To improve the visibility of the features, we dis- 
played the same line-core velocity and line-core intensity maps 
twice in 1st and 2nd column, as the ove rlaid contour lines un- 
avoidably tend to hide the fine- structure. iLeenaarts et ail d2009l) 
found a similar relation between velocity and intensity for the 
Call IR line at 854 nm in observations and simulations. 




Fig. 5. (1st column, top to bottom): line- 
core velocity of Call H, line-core intensity at 
396.849 nm, intensity of H 2 y emission peak. 
(2nd column, bottom): H 2 r. The upper two 
panels are identical to 1st column, but with 
overlaid contours. (3rd column, top to bottom): 
wavelength-integrated polarization degree p toU 
H-index, number of emission peaks in spec- 
trum (black/dark grey/light grey/white corre- 
spond to 0/1/2/3 peaks). Red contours outline 
strong polarization signal, blue contours high 
line-core intensity, and yellow contours large 
blue shifts of the Call H line core. The green 
rectangle/parallelogram denote the locations 
of fibril-like structures. Green contours m Pt t 
mark profiles with two emission peaks. 



The apparent intensity increase in the line-core intensity map 
does not extend far beyond the locations of large photospheric 
polarization signal (blue and red contour lines, respectively), 
whereas for the H-index the brightest parts are to first order co- 
spatial with the polarization signal (middle right panel). The H- 
index is not sensitive to the Doppler shifts producing the bright 
fibrils due to its extended wavelength range, contrary to the line- 
core intensity map. 

3.2. Direct signature of magnetic fields in spectral shape 

Magnetic fields leave, however, also a more direct signature 
in the shape of the core of the Call H line profiles. A spec- 
tral pattern with double-peaked emission on the locations of 
the H2R and H2V emission peaks is found almost exclusively 
together with photospheric magneti c fields dLiu & Smith! Il 972; 
Reza ei et al. 2007; Beck et al. 2008). But even if both emission 
peaks are present at the same time in many cases, they usu- 
ally show a pronounced asymmetry with higher intensities in 
H2V. This asymmetry indicates t he presence of mass flows (e.g., 
lHeaslevl[T975l:lRezaei et al.ll2QQ7h . 

To determine the locations inside the FOV with a direct mag- 
netic field signature, we have analyzed the shape o f the Ca II H 
profiles with a routine similar to the one used by Rezaei et al. 
(2008); the routine searches for local intensity maxima in the 
spectrum. Appendix [A] explains the procedure in more detail 
and shows several examples of its results on different profiles. 
With this routine the intensity of H2V or H2R can only be de- 
termined wh e n the peaks are present in the spectrum. Following 
Reza ei et al. 1 d2QQ7h we thus also determined the integrated in- 
tensity in two 27 pm wide wavelength bands located symmetri- 
cal around the line core that yield an estimate of H2V and H2R 
intensity for each profile. 

The bottom row of Fig. [5] shows the results of both ap- 
proaches that use the shape of the spectrum near the line core. 
The appearance of spectra with both H2V and H2R emission 
(green contour lines in the polarization degree map at upper 



right) is basically restricted to the locations with photospheric 
polarization signal, with an additional extent of about 5" dis- 
tance around the photospheric fields. In the rest of the FOV, 
either single peaked profiles or those without intensity reversal 
are prevailing. The area fraction of double-peaked profiles in the 
FOV is 13%, slightly larger than the estimate of photospheric 
magnetic field coverage. The maps of H2V and H2R using the 27 
pm wavelength bands are displayed in the bottom row, with an 
identical display range; they show that the H2V emission is usu- 
ally stronger than H2R, even on the locations with photospheric 
magnetic fields. For half of the profiles with two emission peaks 
the asym metry between H2V and H2R is above 10 %. As dis- 
cussed by lRezaei etafl (l2007h and lBeck etaD <|2008), this is re- 
lated to the fact that the Ca profiles even above photospheric 
magnetic fields contain a (dominant) contribution with the sig- 
nature of ac oustic shocks, in dicated by a pronounced H2V/H2R 
asymmetry dBeck et alJl2008l their Fig. 17). 

In total, the area fraction with chromo spheric magnetic field 
signature is small (-15 %), half of the affected profiles show in- 
dications for the presence of flows as well, and we are interested 
in the velocity statistics of all type of waves regardless whether 
they are related to magnetic fields or not. We thus assume that 
the magnetic field topology can be neglected in the analysis of 
the velocity statistics of the full FOV. To investigate whether the 
chromospheric magnetic fields have any effect on the LOS ve- 
locities, we will use the statistics of all locations with a double- 
peaked emission in the large-area map (light grey area in bottom 
right of Fig. [5]) later on for a comparison with the full FOV. These 
locations cover most of the photospheric polarization signal as 
well, and thus indicate magnetic fields at all heights. 

We would like to also point out that the bright fibrils may 
nonetheless actually still outline chromospheric magnetic field 
lines. The reason for their appearance may be, however, not the 
presence of the field lines alone, but rather the passage of a wave 
or mass flows along them. This corresponds to one of the indirect 
magnetic heating scenarios discussed in the introduction. 
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Fig. 6. Left column: scatterplot of continuum intensity vs LOS velocities. Dash-dotted lines show a linear regression fit. The data 
has been binned. Right column: scatterplot of the LOS velocities vs their respective line-core intensities. The green pluses show the 
same for locations with two emission peaks in Ca II H. 



4. Intensity-velocity relation 

To obtain the intensity maps (outer wing (OW), middle wing 
(MW) and inner wing (IW), H-index) in the left columns of 
Figs. [3] and HI that should correspond to approximately the same 
formation height as the LOS velocities, we averaged the inten- 
sity in the Ca II H spectrum over the line cores of the respective 
lines (pluses in Fig. [2]). Even if the averaging extends over some 
spectral pixels, we will address the resulting intensity maps as 
the "line-core intensity" of the respective line in the following. 
The approach, however, only partly succeeded in the desired ef- 
fect of producing intensity maps looking similar to the velocity 
maps. For the lowest forming spectral line (line 1, display range 
± 1 kms -1 ; Fig. |4]), the spatial structures in the velocity and the 
line-core or continuum intensity map do match: white in veloc- 
ity (= blue- shifts) corresponds to an increased intensity, but for 
the other lines the relation between intensity and velocity is less 
clear. The display range was increased to ±1.5 kms -1 for line 9 
and ±5 kms -1 for line 8. For all velocity maps besides vl we 
reversed the velocity display range (white = red- shifts). 

In Fig. Owe tried to quantify the relation between intensity 
and velocities for two cases: velocity vs the continuum intensity 
I c (left column) and velocity vs the respective line-core intensity 
of the spectral line (pluses in Fig. [2 right column). The line-core 
intensities were normalized with their mean values. To improve 
the visibi lity of the depende nce, the data were binned in the same 
way as in lBeck et al.l (120071) . We took the pairs of corresponding 
(/, v) points, sorted the intensity values into fifty equidistant bins 
[/, / + A/], and then used the average values 1 = a/] and 

v = (v)[ij+m] in each bin. 

For lines 8 and 9, no correlation with the continuum inten- 
sity remains, whereas line 5 and line 1 show a weak and strong 
anti-correlation, respectively, with the granulation pattern of I c . 



Comparing the velocities to the respective line-core intensities, 
already line 5 shows a reversal of the relation, with red-shifts 
corresponding to increased intensities. Restricting the FOV to 
only the locations with double-peaked emission gives essentially 
the same relationship between velocities and intensities (green 
pluses). 



A second peculiarity of the velocity and intensity maps re- 
lates to the spatial scales. Whereas for the velocities the spatial 
scales of structures in the large-area map (Fig. [4]) increase to 
a more diffuse pattern, this trend is less prominent in the line- 
core intensity maps. Along the temporal axis of the time series 
(Fig. [3}, the velocity patterns do not turn diffuse with height. 
This should be related to the frequency shift of the oscillations 
with height (cf . Appendix |B} that lead to a compression of the 
wave pattern in the temporal domain. 



In brief, both the intensity and the velocity maps other than 
for line 1 do not show any clear relation to the granulation pat- 
tern, implying that they are dominated by processes other than 
the convective motions. The spatial structures in the velocity 
maps show a typical spatial scale that is larger than individual 
granules and that could be p art of the inverse granulation pat- 
tern (e.g. iRutten et al.ll2004l) . These authors concluded that the 
inverse granulation pattern should be a mixture of convection 
reversal and a contribution from gravity waves whose relative 
importance increases with height. The little visual correspon- 
dence between velocity and intensity maps then again supports 
the claim that line 5 and those higher up are dominated processes 
unrelated to convection. 



Fig. 7. Histograms of LOS velocities. Spectral line and width of 
the Gaussian fit (red line) are indicated in each plot. The blue 
curves and values at right refer to the Fourier- filtered velocity 
maps of the time-series; the green curves and the values at left 
are for locations with two emission peaks in Call H. 




Fig. 8. Fourier-filtered velocities of line 1 and 5 of the time- 
series. The grid lines have a spacing of 10" (y-axis) and 7 min 
(x-axis). 



5. Statistics of LOS velocities with height 

The histograms of the line-core velocities in Fig. [7] clearly show 
that the width of the velocity distribution increases with height, 
as the lines have been sorted according to their formation heights 
given in Table [2 We made a Gaussian fit to the distributions to 
define the root-mean- square (rms) value <r Vi of velocity varia- 
tions from the width of the distribution. For the data from the 
time-series, it is possible to remove the velocities induced by 
the convective motions by Fourier filtering. We took out all ve- 
locities corresponding to frequencies below 1.3 mHz (-760 sec 
period) from the velocity maps of line 1, 5, and 9. Figure[8] shows 
the resulting velocity maps for line 1 and 5. The removal of the 
low-frequency variations changes the distributions of velocities 
slightly by reducing their width by 60 (20) ms" 1 for line 1 (line 
5). For line 9, which showed no correlation to the granulation 
pattern, the effect is negligible (blue curves in Fig. [7]). In the 
following, we only used the low-pass filtered velocity from the 
times-series for line 1 and 5, neglecting the large-area map for 
these lines, but used both data sets for lines 8 and 9 where the 
filtering is negligible. Restricting the FOV again to only the lo- 
cations with double-peaked emission (green lines) yields distri- 
butions with a slightly reduced width for all lines. We will not 
use the statistics for the magnetic locations further; the main dif- 




Fig. 9. Relation between co-spatial velocities of lines 1,5,9 and 8 
in the time- series. Top to bottom: scatterplots of vl vs v5, vl vs 
v9, and vl vs v8. Black dots: all data points; red: after binning. 
Blue line: linear regression to the binned data. The slope is given 
in the upper left corner in each scatterplot. 



ference to the full FOV is the reduced rms velocity value for 
Call H, whereas for the other lines the reduction is less severe. 

Up to this point, three arguments show that the patterns in 
the velocity maps cannot be of convective origin: 

- The low-forming lines have been Fourier- filtered for the 
low-frequency part which removes the granular contribution; 
only the filtered time-series is used for lines 1 and 5. 

- The velocity maps for lines 9 and 8 have no correlation to 
the granulation pattern. 

- All lines are dominated by oscillatory behavior with recur- 
rent positive and negative velocities with periods shorter than 
granular time scales (cp. Figs. [3] and [8]). 

The patterns thus should be related to waves of some kind; using 
the velocity maps of the time-series shows that they are propa- 
gating in height as well. 

The question of the wave energy during the propagation 
through the atmosphere is related to the growth of the velocity 
amplitude of the vertical LOS velocities with height that has to 
compensate the reduction of gas density. Figure [9] compares co- 
spatial and co-temporal velocities in the time- series. To obtain 
the velocities co-temporal to line 1, the velocity maps of line 5, 
9, and 8 were shifted in time by 5 sec, 10 sec and 59 sec, re- 
spectively. Figure [TO] shows that with these temporal shifts, the 
velocities of the lines at one location along the slit are well in 
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Fig. 10. Example of LOS velocities at a fixed position along the 
slit in the time-series. Black, blue, green and red lines show the 
velocity of lines 8,9,5, and 1, displaced by +4,0,-2,-3 kms -1 . 
The dashed lines mark the respective zero velocity. The colored 
crosses result from multiplying vl with the coefficients of Fig. [9] 
for v5 and v9, and the ratio of the rms velocities for v8, respec- 
tively. 

Table 2. Line-core formation height (2nd column) and velocity 
amplitude ratios of left to right observed rms velocities, velocity 
of vl propagated in linear wave theory, same for v5, and from 
the fit to the scatterplots. 

line z [km] cr VJ -/cr vl (r v i tProp J(r v i Q- v5fProp J cr v i fit 



1 250 1 1 - 1 

5 450 1.28 1.65 1.28 1.25 

9 600 1.99 2.65 1.93 1.66 

8 1000 7.37 8.75 8.83 1.23 
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Fig. 11. Top: statistical properties of LOS velocities: minimum, 
maximum and rms value (red, black, green). The horizontal dot- 
ted lines are at ± 1.5 kms -1 , the dashed at zero. The horizontal 
error bars correspond to a ±75 km error in the assumed forma- 
tion heights. Bottom: comparison between observed rms veloci- 
ties (green) and those from a propagating of velocities in height 
using linear wave theory. Blue line: propagation of vl; purple 
line: same for v5, red line: same for v9. 



phase for lines 1,5, and 9. For the Ca line core, the agreement of 
phase is worse, but for most of the large- amplitude velocity ex- 
cursions the 1-min shift gives a good agreement (e.g., at t= 5, 3 2, 
or 40 min). For the co-spatial and co-temporal velocity maps, 
we did scatterplots vs the velocity of line 1 (Fig. [9]). The differ- 
ent lines show a strong correlation with the lowermost forming 
line 1 with an ever increasing velocity amplitude ratio. The ratio 
was determined by a linear regression fit to binned values of the 
scatterplot (red dots). The obtained ratio was then used to scale 
up the velocity of line 1 to the velocities to be expected for the 
other lines (colored crosses in Fig. [TO]). For line 5 and line 9, the 
usage of a constant scaling factor reproduces the actual observed 
velocity amplitudes well, regardless of the period of the veloc- 
ity oscillations. For the Ca line core, the scaling coefficient from 
the scatterplot (~1) leads to much too small predicted velocities; 
we thus scaled line 1 up with the ratio of the rms velocities in- 
stead (~7). This reproduces only the range of observed Ca line 
core velocities; no close match in phase or amplitude at a given 
time is achieved. In brief, all velocity oscillations of line 1 ap- 
pear in the other lines with both a time-lag and and an amplitude 
increase. 

The statistical properties of the lines are summarized in the 
upper part of Fig.Qj] the extreme velocities increase from ±1.5 
kms -1 in line 1 to above ± 8 kms" 1 for the line core of Call H, 
the rms value increases from 0.3 kms -1 to 2 kms -1 . 

To investigate the consistency of the observed velocities, we 
propagated the velocities in height using the linear perturbation 
approximation for wave propagation. We used the Fourier trans- 
form of the velocity maps of lines 1, 5 and 9, and propagated the 



velocities upwards using (Mihalas & Weibel Mihalas 1984) 



v(z + Az, cS) = Vi(z, cS) • exp 



1 

2H 



Az 



(2) 



A constant scale height of H = 100 km was used, the sound 
speed c s was set to 7 kms -1 and an acoustic cutoff frequency of 
v c = 5 mHz was used to calculate oj c = 2nv c . For each height 
Az we then calculated the velocity histograms and rms velocities 
like for the observations (see Appendix [B] for more details). The 
lower part of Fig. [TT] compares the observed and synthetic rms 
velocities from the propagation. The propagation of vl leads to 
rms velocities larger than the observed ones for the other lines, 
but if v5 is propagated with Eq. (|2]) a good match to the observed 
values of line 9 and 8 is obtained. This indicates again that the 
velocities from line 5 and higher up are dominated by waves. 
Table [2] compares the observed velocity amplitude ratio between 
line 1 and the other lines with that of the propagated velocities 
and the linear fit to the scatterplots of Fig. [9] Besides the values 
from the propagation of vl, the other values agree well that the 
velocity amplitude increases on average by (1.3, ~2, ~ 8) for 
-=(450, 600, 1000) km. 

Comparison to a semi-empirical atmosphere model All ve- 
locities can be converted to a mechanical energy equivalent by 

.,2 



■ j mech 



= l/2p. 



gas 



(3) 



More precisely, E mec h is the ram pressure and can be considered 
as a mechanical energy density. To distinguish it from the mass 
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Fig. 12. (Top panel): mass energy density, E/p = 1/2 vf, vs 
height. Red, green and black lines correspond to minimal, rms, 
and maximal velocity, respectively. The blue line gives the en- 
hancement of internal energy Api nt required to obtain the orig- 
inal HSRA temperature stratification. The vertical black lines 
indicate the extension of the velocity distributions for each line 
that are shown in the lower four panels. There the blue vertical 
lines denote the energy required to obtain the HSRA model at 
the specific formation height of the line. 



energy density used below, we will refer to the quantity as "vol- 
ume energy density" in the following. 

For the comparison with a semi-empirical solar atmosphere 
model, we chose the temperature enhancement over radiative 
equilibrium (RE) temperature as the quantity that should cor- 
respond to the volume energy density. RE conditions can be 
thought of as providing a lower limit to the temperature for the 
upper photosphere and low chromosphere on a long-term av- 
erage because of the pervasive radiation field from continuum 
layers that c onstantly tries to re-heat plasma above to RE tem- 
perature (e.g. lUlmschneider et alii 197 81: ICheung et al.ll2007|). W e 
selected the HSRA atmosphere model dGingerich et al.l fl97lb . 
because the diffe rences between the various semi-empirical 
models (e.g., in iGingerich et al.l 1 19711: IVernazza et al.l 1 198 it 
iFontenla et aDll999l 120061: lAvrettl l2QQ7h on the location and 
strength of the chromospheric temperature rise are minor in 
the present context. To obtain a RE atmosphere model, we ex- 
changed the temperature stratification of the HSRA atmosphere 
model with that of the Holweger-Mueller model (HOLMUL, 



iHolweger & Mueller] 1 19741) for all optical depths smaller than 
log r=-4 (z~550 km). The HOLMUL model roughly corre- 
sponds to an atmosphere in RE condition (cf. Fig. \T3\ for a plot 
of all stratifications). As the semi-empirical atmosphere models 
like HSRA are assumed to correspond to a temporal and spatial 
average of atmospheric conditions, we will do the comparison to 
the velocity distributions, because their range should cover the 
average status predicted by the semi-empirical model. 

We convert the temperature enhancement over RE condi- 
tions, A r, to an enhancement of the internal energy density, Ap int 
by 

RAT 

kpint = Pgas ' — 7T , (4) 

Mr-i) 

where R = 8.31/mol -1 ^" 1 , ji = \3gmo\~ l and y = 5/3 are 
the gas constant, the specific mass and the adiabatic coefficient, 
respectively. 

For the comparison of Ap int with the mechanical energy of 
the waves, it actually is not necessary to take the gas density p gas 
into account. p gas in chromospheric layers usually comes with a 
large uncertainty: if the density is derived from the condition of 
hydrostatic equilibrium, it strongly depends on the assumed tem- 
perature stratification. One can, however, simply compare 1/2 v 2 
with R/(ji(y- 1)) • AT. Both equations yield a (mass) energy den- 
sity in Jkg -1 . Consider for example the Ca line-core velocity and 
the corresponding enhancement of HSRA over RE at 1000 km. 
R/ji(y - 1) is around 9600 and AT -2000 K; thus, the internal 
energy enhancement is 9600 • 2000 ~ 1.9 • 10 7 Jkg -1 . This is to 
be compared with 0.5 cr 2 8 = 0.5 • 2100 2 ~ 2 • 10 6 Jkg" 1 from the 
LOS velocity, which is an order of magnitude smaller. 

Figure [12] shows how the enhancement of the internal energy 
density required to maintain the original HSRA model above ra- 
diative equilibrium compares with the velocity statistics. For the 
lines 1 and 5, with a formation height below z<500 km, actually 
no additional energy is required to maintain the HSRA model. 
For lines 9 and 8, the extreme velocities (black line and red line) 
would supply more than enough energy compared to the require- 
ments of HSRA, but the rms velocities are clearly insufficient. 
This is highlighted in the distribution of the mass energy density 
in the lower panels. The energy required to obtain the HSRA 
model lies at the upper limit of the distributions for lines 9 and 
8, and is far from the rms value. We point out that the mass en- 
ergy density increases with height, reflecting the increase of the 
velocity amplitude. 

When the gas density is multiplied to obtain a volume energy 
density, the trend with height is reversed. We used a gas density 
with height that corresponds to our RE temperature stratifica- 
tion put to hydrostatic equilibrium; the density is close to that of 
the HSRA model up to z=1000 km. As the histograms and the 
value required to obtain the HSRA would simply be scaled by 
the same coefficient, the same result as in Fig. [El would result 
for the lower panels. We thus have used the plot of volume en- 
ergy density vs height (upper panel of Fig. [13]) differently. We 
made a (manual) fit to the exponentially decreasing volume en- 
ergy density with height to reproduce the observed values (dash- 
dotted lines). For both extreme and rms velocities, a reasonable 
agreement can be achieved with an exponential decay exp(-z/a) 
with a~185 km (middle panel). For easier comparison with other 
work, we converted the volume energy density here to an acous- 
tic flux density by multiplying with a constant sound speed of 
c s ~ 7 kms -1 . We overplotted the curve that comes from propa- 
gating the velocity of line 5 upwards (purple line). The horizon- 
tal dashed line denotes 4 .3 kWm 2 s -1 as a ref erence of chromo- 
spheric radiative losses (IVernazza et al.|[l976h ; for z > 500 km, 
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Fig. 13. Volume energy density, E/V = l/2pv^, vs height (top 
panel). The green and black dash-dotted lines show a fit of an 
exponential decay. {Middle panel): acoustic energy flux in log- 
arithmic scale, obtained as E/V • c s . The purple line shows the 
result of propagating v5 in height. The dashed horizontal line 
denotes a commonly used chromospheric energy requirement. 
The lower panel shows the temperature stratifications resulting 
from the dash-dotted and purple lines of the upper panel in com- 
parison to the original (modified) HSRA and the HOLMUL at- 
mosphere model. 



the energy contained in the LOS velocities falls short of it. If the 
exponential decay curves of the upper panel are converted to the 
corresponding temperature enhancement over RE temperatures 
by equating Eqs. 0]) and ©, the extreme velocities would again 
be more than sufficient for the HSRA requirements, whereas the 
rms velocities would at least suffice for an increase of 500-1500 
K above RE temperatures for heights above 1200 km. 

Thus, in this comparison of mechanical energy contained in 
mass motions and the enhancement of internal energy required to 
lift an RE model to the original HSRA model, the mass motions 
are insufficient to supply the necessary energy. Our interpretation 
of all line-core shifts as being due to mass motions corresponds 
to an upper limit, as some of the shifts of the blends inside the 
Ca II H line may also actually have been due to temperature ef- 
fects. 



6. Summary and discussion 

We derived the mechanical energy contained in the LOS veloc- 
ities of various spectral lines that form in different geometri- 




Fig. 14. Left to right: Ca II H line-core velocity, low-frequency 
contributions, integrated absolute circular polarization signal. 
Velocities are displayed in a ±5 kms -1 range. 



cal heights in the solar atmosphere. The lines could be grouped 
into four height ranges of which only the lowest two show a 
clear relation (anti-correlation) to the granulation pattern in the 
continuum intensity. When comparing the line-core velocities to 
the corresponding line-core intensities, a positive correlation is 
found for all lines but the lowest forming one: red- shifts are re- 
lated to intensity increases. The amplitude of the velocity oscil- 
lations increases with height, both for the whole velocity distri- 
butions and co- spatial (but not co-temporal) velocities in a time- 
series, in good agreement with theoretical predictions from lin- 
ear wave theory. When the wave energy is converted to a cor- 
responding temperature enhancement of the internal energy, it 
falls short of the requirements of lifting a RE atmosphere model 
to the original HSRA model with a chromospheric temperature 
rise. The distributions of the LOS velocities barely reach these 
requirements for the most extreme, and thus rare, velocities. A 
fit of an exponential decay to the mechanical energy density with 
height yields a possible heating and a temperature increase in the 
chromospheric layers of around 500 to 1500 K, which is signifi- 
cantly less than in the HSRA atmosphere model. 



Magnetic field influence on the Ca II H line core The area frac- 
tion of Ca II H profiles that show signatures of magnetic fields 
is around 15 % (Sect. [3}. The rms velocities of all spectral lines 
considered are reduced on these locations, most prominent for 
the Call H line core (Sect. [5]). The line-core velocity of Ca also 
shows a considerable contribution at low temporal frequencies 
(y <1.3 mHz) that is not present in this way for the other lines 
(Appendix |B]). This low-frequency contribution to the LOS ve- 
locity comes m ainly from locations with photospheric magnetic 
fields (Fig.[T4l lLites et al.lll993h . The locations of the magnetic 
fields can nearly be determined from the Ca line-core velocity, 
by selecting the large- amplitude slow variations in Fig. [14] This 
is clear for the two persistent magnetic patches at y=20" and 
50", but also at y=30" a transient and weaker polarization sig- 
nal is seen. Due to the small extent of the spatial scan in the 
time- series we cannot exclude that the polarization signal actu- 
ally comes from a network element that the slit has intersected 
just at the very edge. 

The importance of this field-related velocity component 
should increase with height in the atmosphere when the mag- 
netic field lines spread out to fill the whole chromospheric vol- 
ume. In the formation height of Ca II H, only the locations with 
also photospheric polarization signal are affected which suggests 
that the Ca line in quiet Sun forms below the canopy layer that 



marks the transition from a basically field-free to fully magne- 
tized atmosphere. 

Range of temperature variations If the velocity distributions 
are converted to their (relative) temperature equivalent, the ex- 
treme velocities of the Ca line core lead to a sufficient en- 
hancement of around 2000 K above RE temperature to reach 
the HSRA temperature values. The distributions span the full 
range from to the 2000 K, implying a variation of the temper- 
ature by several hundred K. For line 9 at around 600 km height, 
the temperature variation range al ready covers around 1000 K. 
This is at odds with the claim of lAvrett et al.l (120061) that tem- 
perature variations in the low chromosphere should not exceed 
±200 K away from a "standard" temperature rise comparable 
to that in the HSRA. Their conclusion was based on the small 
amount of intens ity variation in spectra observed with SUMER. 
lAvrett & Loesert (120081) repeat this claim by stating that "the 
brightness temperature variations observed at high spatial res- 
olution, away from active regions, do not exceed a few hundred 
degrees'". They also suggest strongly that any temperature fluctu- 
ations should be around a semi- static model with a temperature 
rise. We will address the question of the range of intensity fluc- 
tuations in the 2nd of the paper of this series in more detail, but 
remark already here that the intensity variations cover a similar 
dynamical range as the velocities and actually imply tempera- 
ture fluctuations above 500 K around a RE atmosphere w ithout 
any chromospheric temperature rise. iRezaei etaLl (120081) have 
shown examples of Ca II H profiles that cannot be reproduced 
with temperature stratifications including a "standard" chromo- 
spheric temperature rise, but actually most of the observed Ca II 
H spectra disagree all the time with the assumption of a strong 
semi- static temperature rise that is only slightly modulated by 
dynamic events. 

Relation between chromospheric intensity and LOS velocities 

The line-core position of Ca II H may be an unsuited tracer of 
the chromospheric dynamics due to the complex line formation. 
We thus use the intensity close to and in the line core of Ca II 
H to discuss the causal relationship between chromospheric in- 
tensity and photospheric velocities. Most of the LOS velocity 
oscillations have temporal frequencies below the commonly as- 
sumed cut-off frequency of v c = 5 mHz (cf. Fig. IB.ll in Appendix 
|B} a nd thu s shoul d not be able to reach chromospheric heights. 
iBeck et aD d2008), however, found that the phase shifts between 
the Ca II H line-core intensity variations and intensity variations 
in the Ca line wing indicate propagating waves down to fre- 
quencies of around 2 mHz, the same as found by Rutten et al. 
(2004) for G-band and Ca line-core intensities. The same ap- 
plies to the phase shifts between Ca line-core inte nsity and pho- 
tosphe ric L OS velocities (ibi d em, th eir Fig. 10 b). ICenteno et al.l 
(12006b and iKhomenko et al] (120081) have shown that the pres- 
ence of radiative energy losses in the atmosphere can lower the 
cutoff frequency. 

That the frequency of the oscillations is of minor importance 
for the wave propagation, is also demonstrated in Fig.[T5j where 
the same spectra and LOS velocities as used for Fig. [9] are dis- 
played, but without applying any temporal shifts to obtain co- 
spatial velocities that are in phase. The position along the slit 
was at y=973; Fig. [3] shows that there was a distance of at least 
2-3" to the nearest location with detected photospheric flux. In 
addition to the line-core velocities of lines 1 (orange line), 5 
(green), 9 (blue), and 8 (black), we plotted also the intensity at 
396.7 nm, around 0.15 nm distance from the Ca line core, and 
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Fig. 15. Top: spectra at one fixed location along the slit with time. 
The short vertical bars at bottom denote lines 1,5,8 and 9, those 
at top 396.7 nm and the H2V emission peak. Middle and bottom: 
velocities and intensities from the spectra above for the time 0-25 
min and 25-50 min, respectively. Black, blue, red and green lines 
show v8, v9, v5 and vl. The thin blue line shows the intensity at 
396.7 nm, the thick purple line that of F^v- Dashed orange tilted 
lines outline the propagation of perturbations with time. 



at the H2V emission peak (upper two lines). The intensity val- 
ues were referenced vs the Ca II H profile that resulted from the 
RE te mperature stratification in a LTE ca lculation with the SIR 
code dRuiz Cobo & del Toro Iniestall 19921) . The line-core veloc- 
ity of Call H, v8, was treated with a running mean over three 
time steps to reduce the scatter and improve the visibility of the 
velocity extrema. 

It can be clearly seen in Fig. [15] that in 90 % of the cases an 
intensity increase of H2V is preceded by the same at 396.7 nm. 
The increase at 396.7 nm can again be traced back in time to a 
velocity oscillation of initially line 1. In some cases, the velocity 
oscillations of lines 1, 5, and 9 appear almost in phase (e.g., at 
t=31 min), but still trigger an increase of intensity at 396.7 nm 
and H2V with a time-lag, indicating a wave propagation with a 
photospheric origin (cf. also Fig. IB. 2b . The inclined dashed lines 
indicate the propagation from a velocity oscillation of line 1 to 
an intensity increase of H2V in around 1 min. Even if most of the 
power is found below 5 mHz, still exactly these waves have a 
signature in the chromosphere in H2V emission. If one would at- 
tribute a dominant period to the wave pattern of line 1 in Fig.[T5l 
it would be closer to 5 min than 3 min (see also Fig. IB. lb . This 
is in some conflict with ICarlsson & Steinl (1 19971) who found in 
their simulations of wave propagation driven by a photospheric 
piston that "runs with only low frequencies (v <4.7 mHz) do 



not produce grains". This could be reconciled with the present 
observations by the fact that the definition of a "grain" comes 
from the common usage of broad-band filter imaging and refers 
to stronger cases of shocks. The long-period oscillations of the 
velocity of line 1 lead to a transient intensity reversal in the Ca 
line core that would most probably not qualify as bright "grain" 
in most cases. 



Energy flux at different heights The energy flux with height 
of Fig. [13] shows a clear decrease throughout the whole height 
rang e. The absolute am ount is very similar to the numbers given 
by IStraus et al.l (120081 their Fig. 3) who, however, applied a se- 
vere filtering in the Fourier space and used only the velocities 
corresponding to gravity waves. Similarly to the latter authors, 
we also find that the mechanical energy at around z = 500 km 
is of the same order as the chromospheric radiative losses. If 
the formation height of line 9 is correctly assumed to be above 
500 km, the mechanical energy seen in the velocity of this line, 
however, already falls below the common chromospheric energy 
requirements. 

What amount of mechanical energy may be hidden from the 
observations ? In the time-series, one could invoke the tempo- 
ral sampling that excludes all oscillations with periods below 
~ 40 sees. The large-area map for comparison is a stochas- 
tic sample of waves at an arbitrary phase. The velocity his- 
tograms of both types of observations were found to be equiv- 
alent, thus, the temporal sampling of 21 sees does not seem to 
have an influence on the observed velocities. This leaves on the 
one hand the spatial resolution, and on the other the extension 
of the height layers that contribute to the respective line-core 
velocity. From the overview maps of the large-area scan it can 
be seen that the velocity patterns in the higher layers are spa- 
tially resolved, because they are more extended than the struc- 
tures in the intensity maps; the typical size of structures also ex- 
ceeds that of the clearly resolved granulation pattern. The typical 
width of intensity contribution functions in the Ca line wing is 
around 100-200 km (cf. Fig.Q]). Assuming that a similar height 
range contributes to each line-core velocity and a sound speed of 
c s = 7 kms -1 leads to the fact that waves with wavelengths be- 
low <200 km, periods below 28 sec, and frequencies larger than 
35 mHz cannot be observed. The power at frequencies above 20 
mHz can be a ssumed to be negligible for chromospheric heat - 
ing (Fig. EQ iFossum & Carlssonl 120051: ICarlsson et al.l l2QQ7h . 
ev en if this was qu estioned by IWedemeyer-Bohm et al.l ( 2007) 
or Kalkofenl (120071) . Note that all these authors did not con- 
sider oscillations below the cutoff frequency, whereas there are 
some indications that even these waves propagate vertically 
dCenteno et aD 120061: iBeck et al.l [2008: Kho menko etaD 12008: 
Stra us et al.l 120081 or the present Fig. [15]) and thus could con- 
tribute significantly to the chromospheric energy balance. 

Another source of energy that could be missing in the ob- 
servations are wave-related processes that have no line-of-sight 
component at disc centre. This will apply mainly to wave modes 
of the magnetic field lines. The field lines of strong photospheric 
flux concentrations are close to vertical, and several wave modes 
like sausage or kink modes have their main component perpen- 
dicular to the field lines. These waves then will only contribute 
to the statistical analysis of the LOS velocities with their small 
vertical component and remain otherwise undetected. 



7. Conclusions 

We have analyzed the mechanical energy contained in line-of- 
sight velocity oscillations of spectral lines forming at different 
heights. The mechanical energy transported by the waves falls 
short of the requirements to lift a radiative equilibrium model 
atmosphere to one of the commonly used semi-empirical tem- 
perature stratifications with a strong chromospheric temperature 
rise. We find, however, that the chromospheric intensity varia- 
tions are coupled to the LOS velocity of the lowermost forming 
spectral line, with a time-lag of around 1 min, and thus are trig- 
gered by photospheric motions. In the second paper of this series 
we will thus investigate, if the energy in the intensity variations, 
and in the intensity profiles as a whole, agree at all with the semi- 
empirical model. 
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Fig. A.l. Examples of the analysis of the number of reversals 
in Ca II H core. Dashed lines mark the analyzed wavelength 
range. Dotted vertical line: Call H core position. Solid vertical: 
locations of intensity reversals. Wavelengths are /i-396 nm. 
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Fig. B.l. Power spectra of the LOS velocity of (left to right) line 
1, 5, 9 and 8. The power is scaled between minimum and max- 
imum in arbitrary units. The dash-dotted line denotes 5 mHz. 



Appendix A: Analysis of Ca II H profile shape 

The Ca II H line profile can show any number between zero or up 
to three intensity reversals near the line core. The number of re- 
versals can be used for a qualitative estimate of which processes 
have led to the emission pattern. A symmetric double-peaked 
shape with equal intensity in blue (//2v) and red (//2r) emission 
peak indicates a static temperature rise in the absence of strong 
mass flows. Mass flows with vertical gradients lead to an asym- 
metry between H2V and H2R by shifting the absorption core. The 
absence of intensity reversals indicates a temperature stratifica- 
tion without a steep temperature rise in the upper atmosphere 
layers. This can either be a case of a monotonical decrease of 
temperature with height, but also happens during the passage of 
a shock wave when it still is located low in the atmosphere. A 
shock front lifts the intensity symmetrically to the line core un- 
til the reversals of #2v and H2R disappear, just before the onset 
of the strong emi ssion phase dLiull 1 9741; iLiu & Skumanich! 1 9741; 
iBeck et al.ll2008h . 

We used a similar routine as Rezaei et a D (120081) to deter- 
mine the number of emission peaks in the line core. The rou- 
tine passes in wavelength through every profile, and searches for 
local maxima whose intensity exceeds that of the neighboring 
wavelength points. Due to the presence of noise, the spectra are 
treated with a running mean over five wavelengths points prior 



to the search for local maxima. Figure [AJ] shows nine randomly 
picked profile examples. The solid vertical lines mark the lo- 
cations of the intensity reversals found in each profile. The up- 
permost right profile is a typical example for the presence of a 
shock front in the lower atmosphere (lifted wings at H2V and 
H2r); the lowermost left one for a shock in the upper atmosphere 
(H2V emission). The lowermost middle pan el shows a case of a 
reversal-free profile with a low H-index (see Reza ei et al . 2008). 



Appendix B: Fourier analysis and vertical 
propagation of velocity distributions 

Fourier power The time-series of Call H spectra allows for a 
Fourier analysis to determine the dominant frequency contribu- 
tions to the variations of the LOS velocities. Figure IB. II shows 
a k-y diagram of the unfiltered velocity maps of line 1,5, 9 and 
8. The contribution at low frequencies (< 1 mHz) is seen to di- 
minish up to line 9, but is the st rongest contributio n for the line 
core of Call H (cf. Fig.[T4land lRutten et al.ll2004l their Fig. 6). 
The power at medium frequencies (v > 2 mHz) is seen to shift 
towards higher frequencies in lines 1,5, and 9. The greatest part 



of the power is located below the commonly assumed cut-off 
frequency of around 5 mHz {dash-dotted horizontal line). 

Phase differences In Fig. IB. 21 we show the phase difference 
between the lowermost forming line 1 and all others. The veloc- 
ity variations show usually only a small or zero phase difference 
between the different lines, which would be expected if their fre- 
quencies are below the cut-off frequency. No clear relation with 
the line-core velocity of Ca II H can be seen at all. The phase- 
relation between the Ca II H line-core velocity and the low pho- 
tospheric velocities shows a much larger scatter than the phase 
relation between the inte nsity of the H 2V e mission peak and pho- 
tospheric velocities (cf. iBeck et al.ll2008l Fig. 10b). This may 
be caused by the complex formation of the Ca line that renders 
the interpretation of the location of minimum line core intensity 
as a velocity doubtful (see for example the lower left panel of 
Fig. IA. lb . and the fact that the main signature of the waves in 
the chromospheric layers is not a Doppler shift, but an intensity 
enhancement. 
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Fig. B.2. Phase differences of LOS velocities as function of os- 
cillation frequency. Left to right: vl to v5, vl to v9 and vl to 
v8. The relative occurence of a phase difference value is given 
by the color code in the range of 5 % to 20 % (7.5 % to 15 % 
for vl-v8). The red curve gives the center of gravity of the phase 
distribution at each frequency; red vertical dashed lines denote 
frequencies of 2 and 5 mHz, respectively. 
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Fig. B.3. Propagation of the velocity of line 5 to the height of 
line 9. Left to right: LOS velocity map of line 5, its Fourier 
transform, Fourier transform at the new height, velocity map 
at the new height, observed LOS velocity map of line 9. 
Velocities/amplitudes are in kms -1 with the color coding at their 
top. 



The result is transformed back to yield the propagated velocity 
map at the new height (4th column). For comparison, the actual 
observed velocity map of line 9 is show in the 5th column. In a 
direct visual comparison the propagated and observed velocity 
map are identical. 

This close agreement also holds for the distribution of the 
velocity power with frequency (cf. Fig. IB.4I) . The propagated 
velocity of line 5 is again nearly indistinguishable from the ob- 
served one of line 9. For line 9 and for line 8 (Call H line core), 
we used here the velocity maps without filtering for granular low 
frequency power. Line 8 shows a strong contribution from low 
frequencies that is not reproduced by the propagated velocities, 
even if the rms velocity values of both observed and propagated 
velocities turn out to be similar. 



Propagation of velocity with height The velocities of the dif- 
ferent lines can be propagated with height using the linear wave 
theory. The corresponding equation for the propagation con- 
tains a frequency-dependent am plitude scaling and is given by 
(Mihal as & Weibel Mihalas1 fT984) 



v(z + Az, oj) = v(z, oj) • exp 



2H c s ) 

for oj < cl) c , and by 



(B.l) 



v(z + Az, oj) = v(z, oj) • exp |^ • zj for oj > oj c . (B.2) 

We used a constant scale height of H - 100 km, a constant sound 
speed c s of 7 kms -1 and an acoustic cutoff frequency of v c = 5 
mHz to calculate oj c = 2nv c . The frequency-dependent increase 
of the amplitude shifts the distributions of Fourier power towards 
higher frequencies. To obtain the velocity map corresponding to, 
e.g., the velocity of line 5 at the height of formation of line 9, one 
needs the contribution of each frequency to the velocity patterns 
in line 5, v(z V 5,oj). To this extent, we took the Fourier transform 
of each position along the slit in the time-series (Fig. IB. 31 2nd 
column) of the velocity map of v5 (1st column). The velocity am- 
plitudes of the Fourier decomposition are then scaled up accord- 
ing to their frequency with Eqs. (IB. II) and (IB. 2b (3rd column). 
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Fig. B.4. Velocity amplitudes of the Fourier transform, averaged 
along the slit. Top to bottom: line 1,5,9, and 8. Black lines shows 
the distribution in the observations, blue and red lines the one 
resulting from propagating line 1 and line 5, respectively. 



